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FOREWORD 


This report summarizes the results of some additional studies 
regarding the beading and spiking phenomena observed in the M551 
Metals Melting experiment conducted in the M512 Materials Processing 
Facility aboard the Skylab Laboratory. This effort was performed by 
personnel in the Fluid Mechanic Applications Group of the Lockheed- 
Huntsville Research & Engineering Center for the NASA -Marshall 
Space Flight Center. 

This work was performed under the direction of Mr. R. M. Poorman, 
S&E-ASTN-MM, who was the principal investigator on the M551 space 
experiment. 
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ABSTRACT 


A study was made regarding the beading and spiking phenomena 
observed in the M551 Metals Melting experiment conducted during the 
Skylab I mission in June 1973. An analysis was made of the beading 
phenomenon based on the Karman vortex shedding theory. The results 
tend to support the hypothesis that beading which occurred in the 
stainless steel and tantalum samples was a Karman vortex street forma- 
tion. 


A dynamic model of cavity oscillation is discussed to explain the 
spiking phenomenon which was observed in the stainless steel and tanta- 
lum samples. Calculations of spiking frequency indicate that the 
intensity of spiking depends primarily on the vapor pressure and sur- 
face tension properties of the material, and is only slightly affected by 
the level of gravitational acceleration. 
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SECTION I. INTRODUCTION 


The M55I Metals Melting experiment was one of the spaceflight 
experiments conducted in the M512 Materials Processing Facility during 
the Skylab I mission in June 1973. The M512 Facility was an experi- 
mental apparatus designed and developed by NASA for conducting 
materials processing experiments in space. The facility consisted of 
a spherical processing chamber approximately 40 cm in diameter. The 
energy source in the M512 unit was an electron beam with a power 
capacity of between 1.0 and 1.6 kilowatts. A vacuum environment could 
be created by venting the chamber directly to outer space. During the 
Skylab I mission. Astronaut Charles Conrad performed the following 
experiments in the M512 facility: (1) the M551 Metals Melting experi- 
ment; (2) the M552 Exothermic Brazing experiment; (3) the M553 
Sphere Forming experiment; (4) the M479 Zero-Gravity Flammability 
experiment; (5) the M554 Composite Casting experiment; and (6) the 
M555 Gallium Arsenide Crystal Growth experiment. 

Lockheed was responsible for performing studies to analyze 
the results of the M551, M552 and M553 experiments, including con- 
ducting ground-based (1-g) tests and KC-135 (low-g) aircraft flight 
tests. The results of these studies and some other directly related 
investigations have been documented [_1 through 23]. 

The present document reports the results of some additional 
studies regarding the bending and spiking phenomena observed in the 
M551 Metals Melting experiment as well as in some ground-based elec- 
tron beam welding experiments. The M551 experiment consisted of a 
rotating metal disk mounted perpendicular to the electron beam (eb) 
heat source with the eb impingement point located 6 cm from the center 
of rotation. The velocity of the beam relative to the impingement point 
was 1.61 cm/sec with a corresponding radial acceleration of 5.09 x 
10" g. The eb was focused to approximately 0.15 cm diameter with 
a beam power of 1.6 kW. The three materials used in the experiment 
were 2219 aluminum, 321 stainless steel and tantalum. The thickness 
of the materials varied with angular position around the disk. During 
the continuous weld portion of each disk, both full and partial penetra- 
tion of the disk was achieved by using a constant power input but a 
varying disk thickness. For eack disk, the continuous weld was fol- 
lowed by a dwell portion. In the dwell portion of the weld, the disk 
was stationary while the electron beam impinged on the thick segment 
of the disk, thus creating a large molten pool. The electron beam was 
then shut off and the pool allowed to solidify [22], 
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section ii. beading phenomenon 

The M551 Metals Melting experiment was conducted by Astronaut 
Charles Conrad during the Skylab I mission in June 1973. During the 
experiment, motion pictures were taken for each of the three metal 
disks. The motion pictures taken along with the experiment specimens 
were then brought back to earth for analysis. Shown in Figures 1, 2 and 
3 are three photographs taken of the 321 stainless steel, 2219 aluminum 
and tantalum specimens after being subjected to electron beam impinge- 
ment in the M551 flight experiment. For comparison purposes, ground- 
based tests were also made of the same specimens under similar oper- 
ating conditions. The photographs taken for the three M551 materials 
during ground tests are shown in Figures 4, 5 and 6. The most striking 
phenomenon appearing in these photographs is the "beading" effect which 
occurs simultaneously with total eb penetration in the stainless steel and 
tantalum specimens. No beading is exhibited in the aluminum sample, 
however. 

It was hypothesized [20 and 22 J that the forming of beads 
after the passage of an electron beam was related to the vortex shedding 
phenomenon. When a cylinder moves through a fluid above a certain 
speed, the cylinder sheds behind its wake a periodic pattern of vortices 
which move alternatively clockwise and counterclockwise. This vortex 
formation is known as Karman vortex street. A sketch showing the 
streamlines of a Karman vortex street is given in the sketch below. 



Karman Vortex Street (Diagrammatic); Stream-lines Drawn in a 
System of Coordinates Moving with the Vortex Street 

In this sketch, \ is the wave length or vortex spacing and h is the 
separation distance between the two rows of vortices. 
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FIGURE 1. STAINLESS STEEL DISK AFTER M551 FLIGHT TESTS 
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FIGURE 3. TANTALUM DISK AFTER M551 FLIGHT TESTS 
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FIGURE 6. TANTALUM DISK AFTER M551 GROUND-BASED TESTS 


The shedding frequency, f, i.e., the number of vortices shed 
per unit time, is dependent upon the density, p, and viscosity, p, of 
the fluid, the diameter of the cylinder, D, and the velocity, V, with 
which the cylinder moves relative to the fluid. The interrelationship 
between these quantities can be correlated in terms of two dimension- 
less parameters known as Strouhal number, S, and Reynolds number, 
R, which are defined as 


f x D 


R = ^D. . 

M 

The relationship between S and R was obtained experimentally 
by Roshko [24j for a fluid in large extent. Roshko's experimental 
results and empirical equations between S and R are shown in Figure 7. 

The wave length of the vortices, X, is given by 



where v is the velocity of the vortex which decreases with the distance 
from the cylinder. The magnitude of v is only of the order of 4 to 5% 
at a distance of 10 to 20 diameters downstream of the cylinder. The 
value of v is even smaller when the flow is confined by side walls. For 
all practical purposes v can be neglected in Eq. (1). 

With the above basic background information, one can compile a 
list of data regarding the M551 space flight and ground-based experiments 
in view of the vortex shedding hypothesis. The compiled data are tabu- 
lated in Table 1. The manner in which these data are obtained is described 
as follows. 

The two temperatures selected for calculation for each of the three 
M551 specimens correspond to the melting temperature and predicted 
maximum temperature [2 0j. The Reynolds number was calculated 
based on the velocity of the electron beam relative to the disk, the focus- 
ing diameter of the electron beam, and the density and viscosity of the 
molten metal. The viscosity of liquid metals depends strongly on tem- 
perature. The following equation was used in computing the viscosity 
of the molten metals 


pi 


E/KT 

= M 0 « 


( 2 ) 
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Strouhal Number (S) 


.220 



FIGURE 7. Si ku'u'HAL NUMBER AS FUNCTION OF THE REYNOLDS NUM- 
BER FOR VORTEX SHEDDING BEHIND A CYLINDER [24] 
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Table 1 

COMPARISON OF THEORETICAL AND EXPERIMENTAL BEAD SPACING 


M55 1 
Materials 

Temperature 

T 

(°K) 

Reynolds 
Numbe r 

R 

Strouhal 

Number 

S 

Shedding 

Frequency 

f 

(Hz) 

Theoretical 

Bead 

Spacing 

X 

(cm) 

Flight 

Bead 

Spacing 

X 

(cm) 

Ground 

Bead 

Spacing 

X 

(cm) 

321 

Stainless 

Steel 

1700 

35 

— 

— 

— 

0.356 

0.418 

2253 

72.5 

0.15 

1.61 

1.0 

2219 

Aluminum 

916 

50 

— 

— 

— 



2113 

171 

0.185 

1.98 

0.81 



Tantalum 

3270 

121 

0.175 

1.88 

0.86 

1.57 

0.94 

3833 

181 

0.188 

2.02 

0.80 





where M 0 is a reference viscosity, E is the activation energy, $ is 
the universal gas constant, and T is the absolute temperature. The 
values of Mo and E are different for different materials. For example, 
for the stainless steel specimen, g Q and E q are estimated to be [25J 


M o = O' 3 ? Centipoise 

E = 9.9 Kcal/gm-mole-°K 


After the Reynolds number was calculated, the Strouhal number 
was derived from Figure 7. The shedding frequency can then be computed 
from the definition of the Strouhal number. The theoretical bead spacing 
(i.e., wave length) is calculated from Eq. (1) by neglecting v. The data 
listed in the last two columns of Table 1 are the average bead spacings 
obtained by direct measurements from Figures 1 through 6. ° 

It is seen in Table 1 that except for the aluminum specimen for 
which no beads were observed, the predicted bead spacing based on 
the vortex shedding theory and the experimentally measured average 
bead Spacing are in good qualitative agreement. Quantitative agree- 
ment was not expected since the geometry and boundary conditions of 
the M55 1 Metals Melting experiment were not quite the same as those 
for a long solid cylinder moving in a fluid of large extent. 

theoretically predicted shedding frequency for each of the 
three specimens is also compared in Table 2 with corresponding ex- 
perimental results for both space flight and ground-based tests. Again 
good qualitative agreement among the data is obtained, except for the 
aluminum specimen. 

theorized that surface tension must play an important role 
in the formation of beads from the vortices of the liquid metal during 
cooling and solidification. The larger the surface tension, the more 
easily the beads can be formed. The magnitude of surface tension (<t! 

or e three M551 materials at their respective melting temperatures 
are as foil own- ° r 


321 stainless steel, 
2219 aluminum, 
Tantalum, 


cr = 1750 dyne/cm 
a = 737 dyne/cm 
o = 2150 dyne/cm 


Anof /iuV 10 ??? that the surface tension of 2219 aluminum is only 30 to 
40% of that of the other two materials. Thus, it is quite possible that 
uring e light and ground tests with the 2219 aluminum specimen, 
vortices were present but no beads were formed due to insufficient 
surface tension force. 
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Table 2 


COMPARISON OF THEORETICAL AND EXPERIMENTAL 
SHEDDING FREQUENCY 


M55 1 
Material 

Theor etical 
Shedding 
Frequency 

(Hz) 

Flight 
Shedding 
Frequency 
. __ (Hz) . 

Ground 

Shedding 

Frequency 

(Hz) 

321 

Stainless 

Steel 

1.61 

4.52 

3.85 

2219 

Aluminum 

1.98 

— 

— 

Tantalum 

2.02 

1.03 

1.71 
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The above analysis is, of course, qualitative. The motion of the 
molten metal in the wake of an electron beam is influenced by many 
factors. In addition to the parameters involved in the definitions of 
Strouhal and Reynolds numbers, the state of the fluid motion is influenced 
y the thickness of the disk, confining side walls (between molten and 
solid metals), free surface and surface tension effects, as well as non- 
uniform temperature distributions. 

In general, the effect of confining side walls is to increase the 
critical Reynolds number for the onset of vortex shedding. According 
to Taneda {/•bj, when the ratio of the diameter of the cylinder to the 
distance between the walls is increased from 0 to 0.5*, the critical 
Reynolds number is increased from 45 to about 80. Shair et al. [27] 
have also conducted experiments showing the effect of confining walls 
on the stability of the wake behind a circular cylinder. Their results 
indicated a more pronounced effect than that reported by Taneda. This 
may also explain the absence of beading in the aluminum specimen, 
since the calculated minimum Reynolds number for this case is only 
about 50. It appears, however, that confining wall effects alone do not 
explain the quantitative differences between the predicted and measured 
values of A and f as tabulated in Tables 1 and 2. 

A recent literature search on the subject of vortex shedding re- 
vealed little information regarding the effects of the finite length of the 
cylinder in the low Reynolds number flow region. At high Reynolds 

?,T b \ r T S (R r > 1000) * a shorter cylinder results in larger vortex spacing 
L Tu N °. ° rmation is found at all regarding nonuniform temperature 

and heat transfer effects on vortex shedding. A detailed analysis of 

these problems would be very interesting but this is beyond the scope 
of the present study. * 


The ratio of the diameter of the electron beam to the width of the 
welding in the M551 experiments is approximately 0.5. 
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SECTION III. SPIKING PHENOMENON 


One of the undesirable features of an electron beam weld is the 
so-called "spiking" phenomenon which is a kind of ripple formation in 
the weld seam during solidification. In the M551 Metals Melting experi- 
ments, spiking phenomena were also observed in the partial penetration 
welds of the stainless steel and tantalum samples (see Figures 1 and 3). 
Weld ripples were also clearly seen in Figure 5 which is a photograph of 
the 2219 aluminum specimen after ground tests. Another prominent 
phenomenon observable in Figure 5 is the apparent disappearance of 
ripples in the weld region behind the molybdenum film. The molybdenum 
film has a dimension of 0.5 x 0.5 x 0.0005 inch. It was placed on each of 
the M551 samples in the weld path at a location of 240 degrees. The pur- 
pose of placing a molybdenum film at each sample was to analyze fluid 
motion during melting by micro -sectioning of the weld. 

This section contains a brief discussion and several simple calcu- 
lations to elucidate the mechanisms of the spiking phenomenon and to 
suggest possible reasons for the occurrence and disappearance of such a 
phenomenon. 

In an attempt to understand how some electron beam weld defects, 
such as spikes, cold shuts and base porosity, were formed, Tong and 
Giedt [29] have performed X-ray analysis of workpieces during actual 
eb welding process. The results not only verified that electron beam 
forms a cavity, but also that the cavity oscillates in size and shape. 

Based on this experimental evidence, the authors of Reference [ 29 ] pro- 
posed a dynamic model for the formation of the cavity and explained how 
this model could account for the formation of the aforementioned defects. 

It was discovered that the characteristic shapes of an electron 
beam cavity bear an amazing similarity to the cavities associated with 
the entry of projectiles into water. When the electron beam strikes the 
surface of the molten pool, the intense heating causes the material to 
vaporize and the vapor pressure exerts a downward force on the melt. 

The vaporization material creates a hypothetical body which behaves 
like a projectile and penetrates the molten metal. This "projectile" 
continues to penetrate until the hydrostatic forces cause a deep closure 
or necking of the cavity. This interrupts the electron beam and sepa- 
rates the original cavity into two cavities, a surface cavity and a sub- 
surface cavity. The lower cavity then ceases to penetrate since the 
driving force (vapor pressure) is removed and this cavity fills with fluid. 
The inertia of the deep closure causes an upward flow of fluid until it 
is overcome by surface vapor pressure and the hypothetical body forms 
and penetrates again. 
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An example of the proposed cavity oscillation sequence for 1100 
aluminum is shown in Figure 8, which was based on X-ray and radio- 
graph observations of melt cavity. 

The above model for the oscillation of the cavity can adequately 
de scribe how some of the welding defects are formed. A spike is formed 
each time the electron beam penetrates to the base of the weld and the 
severity of the spiking depends on the frequency of the oscillation and 
the welding speed as well as the material being welded. 

When the molten metal flows into the cavity, it can easily trap 
bubbles of contaminant gases at the base of the weld and, if the cooling 
rate is sufficiently rapid, the fluid may freeze before the gases can be 
convected to the surface, thus giving rise to porosity in the fusion zone. 
Similarly, if the cooling rate is sufficiently rapid, then the melt^'~as it 
falls into the cavity, may interrupt the beam long enough for the walls 
of the spike to freeze so that the fluid, upon falling into the spike, might 
not bond, thus forming a cold shut or crack. The alternate penetration 
and closing of the cavity are also very efficient in mixing the melt and 
results in a very homogeneous fusion zone. 

Based on the above dynamic model, Tong and Giedt [ 29 ] derived an 
analytical equation which predicts the spiking frequency. This equation 
is (see also References [2lJ and [22j) 



where 

u) = spiking frequency 
g = gravitational acceleration 
h = maximum depth of penetration 
P v = vapor pressure of the molten metal 

p = density of the molten metal 

a = radius of the hypothetical body (can be taken 
approximately as the radius of the electron 
beam focusing) 

The preceding analysis considered welding a horizontal plate from 
above. The configuration utilized for M551 ground test was that of an 
electron beam impinging on a vertical disk from the side at a nine to 
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FIGURE 8. CAVITY OSCILLATION SEQUENCE FOR 1100 
ALUMINUM [29] 
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twelve o'clock position with clockwise rotation of the disk past the im- 
mobile beam. However, it will be shown below that the effects of g on 
u) is actually quite small. Equation (3) can also be utilized to predict 
the trends of spiking formation and the effects of the degree of superheat, 
surface tension and vapor pressure properties on the severity of the 
spiking in various materials. 

Simple calculations have shown that the spiking frequency is pri- 
marily controlled by the magnitude of the vapor pressure and surface 
tension which are strong functions of temperature or degree of super- 
heat. The effect of g on w is relatively small because the potential 
force p g h is usually small compared to P v or 2or/a. 

In a space flight environment, g is of the order 10" 4 x g earth or 
smaller, and is essentially negligible in these calculations. In the limit 
of zero gravity, Eq. (3) becomes 


u = 



g — ► 0 


( 4 ) 


It is noted that the higher the vapor pressure and the lower the 
surface tension, the larger will be the spiking frequency. On the other 

2 0 " 

hand, when the value of — exceeds P v , there will be no oscillation ac- 
cording to either Eq. (3) or Eq. (4) (the solution becomes imaginary). 

The vapor pressure of liquid metals can be calculated from the 
following equation 


P 


v 


e B-A/T 


(5) 


where A and B are empirical constants determined for different mate- 
rials. The approximate values of A and B for the three M551 materials 
and molybdenum are tabulated in Table 3. 


By using the data in Table 3 and the surface tension data from 
Reference [25J, the values of P v and 2cr/a are calculated at different 
degrees of superheat for the three M551 materials and molybdenum. 
The results are tabulated in Table 4. 


R* view of Eq. (4) and the comments following that equation, several 
observations can be made regarding the calculated data in Table 4. 

• At 100% superheat, spikes will form for all the M551 
materials. 
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Table 4 


CALCULATED RESULTS OF P v AND 2c /a FOR THE M351 MATERIALS 


Materials 

Melt ing 
Temp. 

T m 

(°K> 

No Superheat 

T = T m 

20% Superheat 
T= 1.2 T m 

50% Superheat 
T r 1.5 T 

m 

100% Superheat 
T = 2.0 T m 

At Mo's Melting 
Temperature 
T = 2880°K 

P 

V 

2cr/a 

P 

V 

2c/ a 

P v 

2c/a 

P 

V 

2c/a 

P 

V 

2o/a 

2219 

Aluminum 

916 

8 x 10' 4 

1.97 x 10 4 

0.24 

1.80 x 10 4 

81 

1.54 x 10 4 

2.59 x 10 4 

1.11 X 10 4 

1.46 x 10 7 

0.16 x 10 4 

321 

Stainless 

Steel 

1700 

12.1 

4.67 x 10 4 

752 

4.22 x 10 4 

4.61 x 10 4 

3.56 x 10 4 

2.79 x 10 6 

2.48 x 10 4 

1.20 x 10 5 

3.12 x 10 4 

Tantalum 

3270 

0.67 

5.74 x 10 4 

297 

5.03 x 10 4 

1.74 x 10 4 

4.47 x 10 4 

1.06 x 10 6 

3.62 x 10 4 . 

0.20 

6.00 x 10 4 

Molybdenum 

2280 

40 

1 

6.00 x 10 4 

2048 

5.53 x 10 4 

1.20 x 10 5 

4.84 x 10 4 

6.10 x 10 6 

3.70 x 10 4 

40.0 

6.00 x iO 4 

i 


NOTES: = Vapor pressure In dynes/cm 2 

a = Surface tension in cfyne/cm 
a *= Radius of the electron beam 



• At 50% superheat, spikes will form in the stainless 
steel sample. 

• No ripples would form in molybdenum if it were welded at 
2880°K. (This temperature is greater than 50% super- 
heat for stainless steel and 200% superheat for aluminum.) 

• At a given temperature, the vapor pressure of molybdenum 
is much less than that of stainless steel or aluminum. 

• At a given temperature, the surface tension of molybdenum 
is greater than that of stainless steel or aluminum. 

The above observations lead to a possible explanation as to why 
ripples disappear in the weld zone behind the molybdenum as clearly 
indicated in Figure 5. The above reasoning is based on the assumption 
that when the electron beam passes through the molybdenum film, the 
vapor pressure and surface tension of the melt are controlled by molyb- 
denum instead of the parent material. The validity of this assumption 
requires, however, further study. 

The effect of gravitational acceleration on spiking frequency is 
shown in Table 5 for the stainless steel sample at three different tem- 
peratures. It is seen that g has very little effect on the spiking fre- 
quency. 
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Table 5 

CALCULATED SPIKING FREQUENCY FOR STAINLESS STEEL 


g 

T = 2550°K 

T = 2883°K 

T = 3400°K 

(cm/ sec^ ) 

(50% Superheat) 

(70% Superheat) 

(100% Superheat) 

0 

28.9 

84.0 

468 

980 

26.2 

82.0 

467 


w i* 

Units in cycles per second. 









SECTION IV. SUMMARY 

Analyses of both the space flight and ground-based tests of the M551 
Metals Melting experiments indicate the following: 

• Good qualitative agreement was obtained between the 
measured results of beading frequency and spacing 
and the predicted values based on the Kirman 
vortex theory. This tends to support the hypothesis as 
advanced in Reference [zo] that the beading phenomenon 
which occurred in the stainless steel and tantalum samples 
was a Karman vortex street formation. 

i Surface tension and sidewall effects appear to be responsible 
for the absence of headings in the 2219 aluminum sample 
after electron beam impingement. 

• The dynamic model of cavity oscillation as proposed in 
Reference [ 29 ] can adequately explain the splicing phenom- 
enon observed in the partial welds of the stainless steel 
and tantalum specimens. 

• Calculated results of spiking frequencies for the three 
M551 samples indicate that the intensity of spiking is 
primarily determined by the vapor pressure and sur- 
face tension properties, and only slightly affected by 
the level of gravitational acceleration. 

• The disappearance of the spiking phenomenon in the 
weld zone behind the molybdenum film could be caused 
by the changing of vapor pressure and surface tension 
of the melt due to molybdenum contamination. 
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